Charge and Spin Hall Conductivity in Metallic Graphene. 
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Graphene has an unusual low-energy band structure with four chiral bands and half-quantized and 
quantized Hall effects that have recently attracted theoretical and experimental attention. We study 
the Fermi energy and disorder dependence of its spin Hall conductivity a'^^ . In the metallic regime 
we find that vertex corrections enhance the intrinsic spin Hall conductivity and that skew scattering 
can lead to values that exceed the quantized ones expected when the chemical potential is 
inside the spin-orbit induced energy gap. We predict that large spin Hall conductivities will be 
observable in graphene even when the spin-orbit gap does not survive disorder. 



Introduction- The low-energy band structure of graphene 
consists of four chiral bands that realize (2+1)- 
dimensional relativistic field theory models with par- 
ity anomalies. The anomalies imply unusual spectra 
in an external magnetic field and quantized and half- 
quantized Hall effects^-^. Theoretical interest^ in these 
unusual electronic systems has increased^ recently be- 
cause of experimental progress^, including measurements 
of the anticipated half- quantized quantum Hall effect. 
One particularly interesting observation, due to Kane 
and MeleSii, is that because of a gap produced by spin- 
orbit interactions, the spin Hall conductivity asH of un- 
doped graphene is quantized in the absence of a magnetic 
field. This suggestion is related to recent work on the 
anomalous Hall effect in ferromagnetic metals^ and on 
its paramagnetic cousin, the spin Hall effectSi, in which 
it was suggested that these transport coefficients can be 
dominated by an intrinsic momentum-space Berry phase 
contribution that reduces to quantized values when the 
Fermi level is in a gap. Here we examine how the quan- 
tized spin Hall effect is altered when the Fermi energy 
in the graphene plane is gated into the metallic regime. 
We find that the intrinsic spin Hall effect is no longer 
quantized, that it is enhanced by disorder vertex correc- 
tions, and that in the metallic regime skew scattering can 
potentially lead to parametrically larger spin Hall con- 
ductivities. Because the Bloch state disorder broadening 
in current samples is (according to our estimates) much 
larger than the clean system spin-orbit gap, these results 
are necessary for the interpretation of experiment. Spin- 
Hall effects should be observable even when the spin-orbit 
gap does not survive disorder. 

Disordered Graphene Model- When spin-orbit interac- 
tions are included,^ the low-energy physics of a clean 
undoped graphene crystal is described by an eight-band 
envelope function Hamiltonian 



^^0 = v{KT^ax + kyay) + Act^t^s^ 



(1) 



where = ± is the up/down electron spin component 
perpendicular to the graphene plane, = ± is a valley 
label that specifies one of the two inequivalent {K and 
K') points in the crystal Brillouin zone near which low- 
energy states occur, and the ai are Pauli matrices rep- 
resenting a pseudo-spin degree of freedom corresponding 



to the two sites per primitive cell of a hexagonal lattice. 
The parameter A is the strength of the spin-orbit cou- 
pling and we take h = 1. For A = this Hamiltonian 
defines four spin-degenerate gapless bands in which the 
pscudospin orientation lies in the x-y plane and winds 
around the z-axis, either clockwise or counter-clockwise, 
with a 27r planar wavevector rotation. The operators ct^, 
Sz and Tz commute with each other. Random defects 
can in general produce transitions between bands and 
between spins. Here we assume spatially smooth spin- 
independent disorder so that Sz and Tz are good quantum 
numbers, allowing us to consider the cases rz,Sz = ±1 
independently. For this disorder model we evaluate the 
Kubo-formula Hall conductivity in the self-consistent 
Born approximation (SCBA) for chemical potentials in- 
side and outside the spin-orbit gap, including both non- 
trivial pseudospin dependent disorder self-energies and 
ladder diagram vertex corrections. When the chemi- 
cal potential lies in the gap, an elementary calculation 
shows that in the absence of disorder the single-band 
bulk partial Hall conductivity is given exactly by the 
half-quantized Berry phase contribution, ^'8. — (sze^/2/i). 
disorder corrections to intrinsic Hall effect are small near 
the gap edge but yield substantial enhancement in more 
strongly gated systems. 

2D-Dirac-band Hall efJect-The 2D Dirac Hamiltonian in 
the spin "f i^T- valley is 



H = v{kx(7x 



kyOy) 



Aa, 



(2) 



Spin-orbit-coupling opens up a gap which breaks the 
spectrum into an electron band at positive energies and 
a hole band at negative energies = i-^/A^ -f {vkY, 
where fc = |k| and ± refer to electron and hole bands re- 
spectively. (The three other graphene bands differ either 
in the Dirac band chirality sense, or in the sign of the 
mass term, or in both ways.) In what follows, we assume 
that the Fermi energy is positive; because of the symme- 
try of the Dirac Hamiltonian generalization to negative 
is trivial. 

The Kubo formula for the Hall conductivity depends 
on both band-diagonal and off-diagonal matrix elements 
of the velocity operator and on the electronic Green's 
function. The disorder-free retarded Green's function 
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and velocity operators for this Hamiltonian are G^(e) = 
(e — H +ir])~^ , = va^, and Vy = VGy. It will prove con- 
venient to use the Streda-Smrckaii version of the Kubo 
formula which separates Fermi surface and occupied state 
contributions: cr^j, = al.y + cr^^ where 



2 f + ca 



47r de 



df{e) 



(3) 



-«,G«(eK(G^(e)-G^(e))] 



and 



7^^ 



2 /• + 00 



47r 



defie)TT[v,G''ie)vy 



de ^ de 



xvyG"ie) ~ v^G^{e)vy^^ + v,^^VyG%m 

2D-Dirac hand intrinsic Hall conductivity — The Hall 
conductivity in the absence of disorder is most simply 
evaluated by expressing^ it in terms of matrix elements 
of velocity operator between unperturbed Bloch states: 



-int e_ 

"xy — n Z^k (f+ 



:2lm[{Ui^\Vy\u+){u^\Va:\Ui^)] (5) 



where the are occupation numbers in the electron 
and hole bands, H. the area of the system, and \u^) the k- 
dependent pseudospinors of the chiral Dirac Hamiltonian, 
Eq.dH. 



cos{9/2) 
sin{9/2)e"^ 



sin{0/2) 
~cos{e/2)e">' 



(6) 

where cos{9) = A/y^wfcp + A^, and tan((/)) = ky/kx- 
For the chemical potential in the upper band with Fermi 
momentum kp we find 



'xy 



(7) 



Eq.(jZj) includes a^y and the disorder free limit of ct^j^. In 
the metallic regime the disorder-independent part of cr^^ 
equals with so that in this regime cr^y = 0. 

When the chemical potential is in the gap (fc_F ^ 0) 



'xy 



xy 



47r 



(8) 



This is the 2D-Dirac model's half quantized (in units 
12'Kh) Hall conductivity, which after summing over 
bands is responsible for the quantum spin-Hall-effect dis- 
cussed in Refs. It should seem surprising that the 
Hall conductivity is only half-integer given general 
arguments that the Hall conductance of non-interacting 
electrons must be quantized. The resolution of this para- 
dox is that bands come in pairs. The sum of the K and K' 
valley bulk conductivities is quantized; correspondingly 
only one band of edge states is induced by the truncation 
of both K and K' bulk bands. 



FIG. 1: Self-energy Feynman diagram in self-consistent Born 
approximation. 



Influence of Disorder on axy '■ — We assume a (5-correlated 
spin-independent random potential with Gaussian corre- 
lations iy{Yi)V{v2))dis = nV^Sin - ra). The SCBA 
that we employ includes only contributions from Feyn- 
man diagrams without crossed disorder correlation lines. 
This common approximation is self-consistent but is in- 
complete. We assume that crossed-disorder-line contri- 
butions give rise to parametrically distinguishable effects 
and do not affect our qualitative conclusions about Hall 
effects in metallic graphene. Fig Q illustrates the SCBA 
self-energy diagram which can be evaluated to obtain 
= — + azC0s{9)) where r'^ is a quantum life 
time at the Fermi surface: 



1/t« =nVjf / kdk5{eF~e+) 



nVi^kp 



(9) 



Following the notation of Dugaev et al}^', the SCBA re- 
tarded Green's function is 



rR — 1 

^ ~ 1/G«-S« 



tF+iVo+v(k^iy^+kyny) + (A—iV-i)G^ 




-4-47") 

(10) 

where Fq ^ 1/(4t9), Fi = Focos(6'), 7+ = Fq -hFi cos(6l), 
7~ = Fq — Fi cos(6'). For these chiral bands disorder not 
only gives the quasiparticle states a finite lifetime but also 
changes the quasiparticle eigenspinors. The SCBA for 



• + 



FIG. 2: Vertex correction Feynman diagram. Black dots rep- 
resent the Pauli operator. 

two-particle Green's functions like the Hall conductivity 
includes in addition ladder diagram vertex correction il- 
lustrated in Fig[5] For large vkpT the terms in axy which 
are products of retarded and advanced Green's functions 
dominate so that the 2D matrix vertex function for which 
we must solve satisfies: 



(11) 



This equation is most easily solved by assuming that 



cay 



da. 



(12) 



and deriving equations for a, b, c, and d. We find that 



_ 2{(vky+2A^) 



8roA(frfc)"-H2A") /-,o^ 
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and that a = d = 0. The SCBA al.y is obtained by sub- 



stituting the disorder-dressed Green's function fEg. ltTUI) ') 
for the bare Green's function e 
formula Eq.@. We find that 



for the bare Green's function and vTy for Vy in the Kubo 



^''y ~ 47r^(vk^y+A^^ ^ iST^JITF ~^ (4A- + («/cF)-)-i- 

(14) 

The second and third terms in square brackets in Ea.f ll4|) 
represent disorder corrections to the intrinsic Hall con- 
ductivity of the 2D-Dirac model. We note that all terms 
are independent of the disorder potential strength and 
of the concentration of scatterers and in this sense are 
parametrically similar. They do however have different 
dependences on the position of the Fermi level. Note 
that when the chemical potential approaches the gap the 
intrinsic contribution remains finite and disorder correc- 
tions vanish, recovering the model's half quantized Hall 
effect. 

Non- Gaussian disorder. — We have so far made the usual 
approximation of assuming Gaussian disorder correla- 
tions. Although normally small, non-zero third mo- 
ments of the disorder potential distribution, can-"^^'^^ al- 
ter axy qualitatively since they can favor scattering with 
a particular chirality (skew scattering) and consequently 
lead to a axy contribution that diverges in the limit of 
weak disorder scattering. The size of this contribution 
to axy is particularly difficult to estimate since it de- 
pends very strongly on the details of the scattering po- 
tential. To illustrate its potential role we consider for 
concreteness a model of uncorrelated (5-function scatter- 
ers: V{r) = ^,F,5(r-R0, R^ random, {V^ = 0, 
{{V^f) ^Vi^Q and {{V,f) = ^ 0. 

Given asymmetric scattering, skew scattering is a more 
physically transparent contribution to the Hall conduc- 
tivity and can be described directly using either Boltz- 
mann transport theory or the Kubo formula, including 
the non-standard Feynman diagrams implied by non- 
Gaussian disorder models. We apply results which have 
been derived previously to the graphene case. Let ijj^ — 
(l/-\/f2)e''"'|uj^) be a Bloch state in the electron band 

with positive energy and Vk,k' = (V'k l^lV'k') be a dis- 
order potential matrix elements within the band. Then, 
following Eqs.(32)-(36) in Ref.l^ for the case of zero 
temperature and a single band we find that 



(er*^)2 



(2^) 



de 



(15) 



where Vx(ii) = de^/dkx, vp is the Fermi velocity, and 

l/r*'^ = /|^^k.k'(l-CO5(0-0')) 



(16) 



Since the scattering rate Wk,k' is usually only weakly 
chiral (r*'' ^ t-^) ti^k.k' can be estimated from time- 
dependent perturbation theoryi^ii^. The lowest order 



symmetric scattering rate is given by the Golden rule ex- 
pression, while the lowest order antisymmetric contribu- 
tion appears at third order (see, for example, Eqs. (2.7) 
and (3.11) in ReflH. 



4S = -i2nme^-e^,) J ^ 

Ini(Vk,k'14;',k"Vk",k)<i-is S{ek. - Ek")- (17) 



This yields 



1 



{vkf 



4A2 



Aivkp)' 



(T«)2nVo^8[(w/fcF)2-f A2]3/2' 



(18) 
(19) 



so that the skew scattering Hall conductivity contribu- 
tion due to non-Gaussian disorder correlations is 



' 27m Vq" (4A2-|-(«fc^.)2)2 



(20) 



The Hall conductivity contribution (|20|l is inversely pro- 
portional to the impurity concentration n, and therefore 
can in principle dominate in the clean limit. Since the 
size of third disorder correlation moment in a particular 
sample is unlikely to be reliably known and can be ex- 
ceedingly small, we expect that the relative importance 
of skew scattering will always have to be assessed exper- 
imentally. 

Application to Graphene- A finite charge Hall con- 
ductance requires broken time reversal symmetry. In 
graphene the vanishing conductance results from cancel- 
lation between bands of opposite spin. The Hall conduc- 
tance we evaluate here could be measured in graphene 
if the Fermi levels in the two spin-]" and the two spin-| 
bands differed. It may be possible to generate spin polar- 
ization in graphene by optical orientation, by tunneling 
through ferromagnetic contacts, or by hyperfine coupling 
to polarized nuclei. We note that the z-component of 
spin is expected to relax particularly slowly in graphene 
because of the planar character of the crystal and the tt- 
character of the orbitals near the Fermi energy. The al- 
ternative of studying the physics we address here, by ap- 
plying an external magnetic field, is not favorable since it 
leads to an ordinary Hall effect in addition to the anoma- 
lous Hall effect. We estimate that the anomalous portion 
of the Hall conductance in an external field is smaller by 
a factor ^ (A/hvk) x (l/vkpr)'^ . When the chemical po- 
tentials of spin-up and spin-down electrons are different 
our Hall effect calculation for each band remains valid. 
The total Hall current is therefore 



a^y""" ^2{axy{fl^)-<Jxy{fil)) 



(21) 



where the coefficient 2 reflects equal contributions from 
the K and K' valleys. 

The Hall conductivity we evaluate appears in the spin- 
Hall response even in the absence of external magnetic 
fields. To find the magnitude of the SHE one should 
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remember that instead of charge e we are interested in 
spin ±1/2 carried by electrons: a^y = Aaxy/2e. Here 
the coefficient 4 is due to the 4 bands which contribute 
equally to the SHE. The spin-Hall effect could be mea- 
sured by using ferromagnetic leads, in the extreme case 
measuring transport only in one spin subsystem. For 
that case the charge Hall conductivity becomes 2axy 
We expect that the results we derive here are valid foB?P 
ep ^ whereas the quantized spin Hall conductivity 
will be observable only if A > r"-'^. The value of r"-'^ in 
current samples can be estimated roughly from measured 
mobilities^ which are roughly constant except for Fermi 
energies below ~ 50meV. Associating the change in mo- 
bility at low carrier densities with disorder mixing be- 
tween electron and hole bands implies a value of the 
same order. The value of A is difficult to estimate accu- 
rately. Based on the relevant potential energy and length 
scales Kane and Mele have estimated that A ~ 0.2 meV. 
This is likely to be an overestimate since the splitting 



represents an average of spin-orbit interactions that vary 
in sign over the system. We2i have separately estimated 
on the basis of a tight-binding model with atomic spin- 
orbit interactions and ab initio electronic structure cal- 
culations that A ~ 0.001 meV. In any event, it appears 
clear that sample quality will need to improve substan- 
tially in order to realize the quantum spin Hall effect. 
As our calculation shows, however, the surprisingly large 
anomalous Hall conductivities that flow from the chiral 
graphene bands should still be measurable in the metallic 
regime. Skew-scattering contributions, if present, should 
be separable experimentally in gated samples on the basis 
of their distinct carrier density dependence. 

Acknowledgments. The authors are grateful for use- 
ful discussions with L. Brey, P. Bruno, F.D.M. Hal- 
dane, C. Kane, Q. Niu, K. Nomura, D. Novikov, and 
S. Urazhdin. This work was supported by the Welch 
Foundation, by DOE grant DE-FG03-02ER45958, and 
by onr-n00014061 0122. 



1 Gordon W. SemenofT, Phys. Rev. Lett. 53, 2449 (1984). 

2 F.D.M. Haldane, Phys. Rev. Lett. 61, 2015 (1988). 

3 D.P DiVincenzo and E.J. Mele, Phys. Rev. B 29, 1685 
(1984); J. Gonzalez, F. Guinea, and M.A.H. Vozmedi- 
ano, Phys.Rev.Lett. 77, 3589 (1996); Yisong Zheng and 
Tsuneya Ando, Phys. Rev. B 65, 255420 (2002); D. V. 
Khveshchenko Phys. Rev. B 65, 235111 (2002); K. B. Efe- 
tov and V. R. Kogan, Phys. Rev. B 68, 245313 (2003); 
Paul E. Lammert and Vincent H. Crespi, Phys. Rev. B 69, 
035406 (2004). 

* V. P. Gusynin and S. G. Sharapov, Phys. Rev. Lett. 95, 
146801 (2005); V.P Gusynin, S.G. Sharapov, preprint 
cond-mat/0512157 (2005); N. M. R. Peres, F. Guinea, and 
A. H. Castro Neto, Phys. Rev. B 72, 174406 (2005). 

^ Novoselov, K. S. et ai, Science 306, 666 (2004); Y. Zhang, 
J.P. Small, M.E.S. Amori, and P. Kim, Phys. Rev. Lett. 
94, 176803 (2005); Yuanbo Zhang, Yan-Wen Tan, Horst 
L. Stormer, Philip Kim, preprint cond-mat/0509355 K.S. 
Novoselov, A.K. Geim, S.V. Morozov, D. Jiang, M.I. Kat- 
snelson, LV. Grigorieva, S.V. Dubonos, A. A. Firsov, Na- 
ture 438 (2005) 197, preprint cond-mat/0509330 
C. L. Kane and E. J. Mele, Phys. Rev. Lett. 95, 226801 
(2005). 

C. L. Kane and E. J. Mele, Phys. Rev. Lett. 95, 146802 
(2005). 

* G. Sundaram and Q. Niu, Phys. Rev. B 59, 14915 (1999); 
T. Jungwirth, Qian Niu, A. H. MacDonald, Phys. Rev. 
Lett. 88, 207208 (2002); D. Culcer, A.H. MacDonald and 
Q. Niu, Phys. Rev. B 68, 045327 (2003); M. Onoda and 
N. Nagaosa, J. Phys. Soc. Jpn. 71, 19 (2002); Z. Fang, N. 
Nagaosa, K.S Takahashi et ai, Science 302, 5642 (2003); 
Y. Yao et al, Phys Rev Lett 92, 037204 (2004). 

® Shuichi Murakami, Naoto Nagaosa, Shou-Gheng Zhang, 
Science 301, 1348 (2003); J. Sinova, D. Gulcer, Q. Niu, N. 



A. Sinitsyn, T. Jungwirth, and A. H. MacDonald, Phys. 
Rev. Lett. 92, 126603 (2004). 
^° L. Sheng, D. N. Sheng, G. S. Ting and F. D. M. Haldane, 
Phys. Rev. Lett. 95, 136602 (2005). 

P Streda, J. Phys. G: Solid State Phys. 15 No 22 (10 Au- 
gust 1982) L717-L721; L Smrcka and P Streda, J. Phys. G: 
Solid State Phys. 10 No 12 (28 June 1977) 2153-2161. 
A. Crepieux and P. Bruno, Phys. Rev. B 64, 014416 
(2001). 

" V. K. Dugaev, P. Bruno, M. Taillefumier, B. Canals and 

C. Lacroix, Phys. Rev. B 71, 224423 (2005). 
" J. M. Luttinger, Phys. Rev. 112, 739-751 (1958). 

J. Smit, Physica (Amsterdam) 21, 877 

(1955); ibid. Phys. Rev. B 8, 2349 (1973). 
^® P. Leroux-Hugon and A. Ghazali, J. Phys. C: Solid State 

Phys. 5, 1072 (1972). 

J. Schliemann and D. Loss, Phys. Rev. B 68, 165311 
(2003). 

^® L. Berger Phys. Rev. B 2, 4559-4566 (1970). 

N. A. Sinitsyn, Qian Niu, Jairo Sinova and K. Nomura, 
Phys. Rev. B 72, 045346 (2005); N. A. Sinitsyn, Q. Niu 
and A. H. MacDonald, Phys. Rev. B 73, 075318 (2006). 
This follows from the observation that the energy separa- 
tion between states of upper and lower bands with the same 
k are mainly of an order of eF and hence for epr >> 1 in- 
terband matrix elements of velocity operator are not renor- 
malized. This may not be true in other applications, for ex- 
ample in the Rashba coupled 2D electron system—, where 
interband separation is of the size of the spin-orbit cou- 
pling and analogous calculations would be valid only in 
much stronger limit Aso » 1/t. 

Hongki Min, J. E. Hill, N.A. Sinitsyn J. Sinova and A.H. 
MacDonald, unpublished. 



